function f = cara_function(v,alpha,c_with,c_without)

%primatives
wealth = 7e4;

%utility with insurance
u_with = mean(-exp(-alpha*(wealth-v-c_with)));

%utility without insurance
u_without = mean(-exp(-alpha*(wealth-c_without)));

f = (u_without - u_with)^2;